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INTRODUCTION 


Accurate  numerical  analyses  of  stress  and  strain  in  steel  pressure  vessels,  including  gun 
tubes,  require  high  quality  uniaxial  stress-strain  data  and  an  accurate  numerical  fit  to  such  data. 

A  detailed  overview  is  given  in  Reference  1.  These  data  then  provide  the  equivalent  stress  input 
to  yield  criteria,  such  as  von  Mises  and  Tresca,  for  the  numerical  solution  of  two-dimensional  and 
three-dimensional  configurations.  A  serious  omission  from  published  numerical  fitting 
procedures  is  any  consistent  and  reliable  method  for  representing  the  behavior  of  materials  that 
exhibit  Bauschinger  effect,  which  in  turn  serves  to  reduce  the  yield  strength  in  compression  as  a 
result  of  prior  tensile  plastic  overload  (ref  2). 

Three  types  of  steel  are  addressed  herein: 

•  Two  A723  steels  with  nominal  (0. 1  %  offset)  yield  strengths  of  1 1 30  and  1 330  MPa. 

•  H Y 1 80  with  nominal  (0. 1  %  offset)  yield  strength  of  1 1 80  MPa. 

•  PH  13-8  Mo  and  PH  13-8  Mo  super-tough  with  nominal  (0.1%  offset)  yield  strengths  of 
1380  and  1355  MPa,  respectively. 

Details  of  the  experimental  testing  procedure  are  given  in  Reference  3. 

The  principal  features  of  a  uniaxial  test  for  gun  steels  are  illustrated  in  Figure  1. 


Figure  1.  Typical  uniaxial  stress-strain  behavior  showing  strain  hardening, 
reduced  elastic  modulus  in  compression,  and  Bauschinger  effect. 
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LI  represents  an  initial  tensile  loading  regime,  0-A,  during  which  the  steel  behaves 
elastically  up  to  the  yield  point  ay  defined  by  a  given  percentage  offset.  The  elastic  modulus  over 
this  range  is  The  material  then  behaves  plastically,  A-C1-C2.  This  phase  may  involve 
significant  nonlinearity,  A-Cl,  and  predominantly  linear  strain  hardening  (or  softening),  C1-C2. 

In  the  case  of  A723  steel  regime,  A-C2  is  characterized  by  modest  linear  strain  hardening  over  a 
strain  range  of  In  the  case  of  HY180,  it  is  characterized  by  nonlinearity  followed  by  zero 
strain  hardening,  reaching  a  constant  value  at  approximately  3%  total  strain,  whereas  13-8  steels 
exhibit  nonlinearity  followed  by  strain  softening. 

U1  represents  reversal  of  loading  with  elastic  modulus  this  elastic  unloading  with 
modulus  up  to  17%  lower  than  is  discussed  in  Reference  3  and  is  of  importance  because 
it  will  contribute  significantly  in  subsequent  autofrettage  residual  stress  calculations.  The  tensile 
(positive)  stress  reduces  and  subsequently  becomes  compressive  (negative)  at  point  D  but 
continues  to  behave  elastically  up  to  point  E.  Thereafter  behavior  becomes  nonlinear,  moving 
asymptotically  towards  a  bound.  The  value  of  stress  at  point  E,  the  onset  of  nonlinearity,  is  -fiGy 
where  P  is  the  Bauschinger  effect  factor  (BEF)  (ref  2).  Both  P  and  the  shape  of  the  curve  E-F  are 
a  strong  function  of  the  maximum  initial  plastic  strain 

APPLICATION  TO  GUN  TUBES 

When  a  tube  undergoes  autofrettage,  the  equivalent  stress  within  the  gun  tube  follows  0-A- 
C1-C2  during  the  initied  autofrettage  pressurization  or  swage,  and  C2-D-E-F  during  the  removal 
of  swage  or  pressure.  Hence  a  family  of  uniaxial  cycles,  0-A-C1-C2-D-E-F,  each  a  function  of 
initial  plastic  strain  and  hence  of  radial  location,  defines  equivalent  stress  for  the  gun  steel  during 
the  autofrettage  process.  This,  in  turn,  with  appropriate  equilibrium,  compatibility,  and  boundary 
conditions,  is  sufficient  to  calculate  numerically  the  residual  stress  locked  into  the  tube  by 
autofrettage.  Results  of  such  calculations  are  well  documented  for  the  current  range  of  gun  steels 
(refs  4,5),  based  upon  A723  1100  MPa  steel,  and  indeed  conform  extremely  well  to  the  ASME 
PVP  code  (ref  6).  It  is,  however,  necessary  to  determine  whether  the  model  fits  current  gun  steels 
in  all  respects  and  whether  other  candidate  steels  exhibit  similar  uniaxial  behavior.  These  two 
issues  are  examined  within  this  report  and  an  improved  model  is  proposed.  Subsequent  work 
will  address  the  numerical  modeling  and  quantification  of  residual  stresses  in  autofrettaged 
pressure  vessels  made  from  these  candidate  high-strength  steels. 

STRAIN  RANGES  IN  TYPICAL  GUN  TUBES  DURING  AUTOFRETTAGE 

Plastic  strains  during  autofrettage  are  a  strong  function  of  the  ratio  autofrettage  radiusAjore 
radius  and  a  weak  function  of  the  ratio  outer  radius/bore  radius  (ref  4).  Consider  a  typical 
hydraulically  autofrettaged  gun  tube  having  outer  radius/bore  radius  =  2  with  70%  overstrain 
(i.e.,  autofrettage  radius/bore  radius  =  1.7).  In  this  case,  using  program  HENCKY  (ref  4),  and 
assuming  typical  properties  for  an  A723-1 100  MPa  steel,  equivalent  strains  at  the  bore  and  at 
points  10%  into  the  wall  thickness  may  be  calculated.  These  are  listed  in  Table  1. 
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Table  1.  Percentage  Equivalent  Strain  Values  During  Autofrettage  Loading  (LI), 
and  Autofrettage  Unloading  (Ul)  at  Two  Locations  Within  the  Tube  Wall 
(Diameter  Ratio  =  2.0, 70  %  Overstrain) 


At  Bore 

10%  Into  Wall 

Elastic  Strain  (0-A) 

0.531% 

0.531% 

Plastic  Strain  (A-C2) 

1.054% 

0.761% 

Elastic  Strain  (C2-E) 

0.621% 

0.628% 

Plastic  Strain  (E-F) 

0.362% 

0.178% 

Note  the  modest,  but  significant,  value  of  unloading  plastic  strain,  approximately  one-third  or 
less  of  the  plastic  strain  during  loading.  This  serves  to  emphasize  the  fact  that  uniaxial  testing, 
which  is  intended  to  be  truly  analogous  to  material  behavior  at  the  bore  of  a  gun  tube,  requires 
only  small  reverse  strains.  This  would  involve  unloading  to  a  point  roughly  one-eighth  of  the 
way  between  points  E  and  F  in  Figure  1.  A  good  engineering  rule  of  thumb,  which  applies  over 
the  maximum  region  likely  to  be  affected  by  further  yielding  during  proof  testing,  is  seen  to  be  a 
total  strain  reversal  (elastic  +  plastic)  equal  to  initial  plastic  strain.  Table  1  thus  provides  uniaxial 
testing  strains  that  model  autofrettage  loading  and  unloading.  Any  subsequent  loading  cycle,  L2 
(e.g.,  proof  testing)  must  obviously  be  initiated  from  a  realistically  located  point  F. 

The  values  and  locations  of  plastic  strain  during  autofrettage  are  illustrated  in  Figure  2  for  a 
typical  tube  with  70%  overstrain.  The  maximum  extent  of  reyielding  during  autofrettage 
unloading  is  also  shown  (24%  of  wall  thickness),  as  is  the  maximum  extent  of  reyielding  during 
typical  proof  testing  (13%  of  wall  thickness). 


Figure  2.  Percentage  initial  plastic  strain  during  autofrettage  loading 
(continuous  line)  and  maximum  extent  of  initial  yielding 
and  reyielding  (aiTows)(diameter  ratio  =  2.0,  70%  overstrain). 
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MODELING  OF  THE  INITIAL  LOADING  PHASE 


Linear  elastic  behavior  in  Figure  1  is  represented  as  O-A.  The  plastic  phase  A-C1-C2  is 
herein  represented  as 


CTpt^/ tTy  =  <1  Tcinh(c  )  +  d  Cpt^  ( 1 ) 

The  first  term  models  the  predominantly  nonlinear  phase  A-Cl,  while  the  second  term 
provides  linear  strain  hardening  or  softening.  Hence,  for  bilinear  materials  with  zero  strain 
hardening  a  =  d  =  Q,fox  bi-linear  with  strain  hardening  a  =  0,  and  for  nonlinear  with  zero  strain 
hardening  d=0. 

MODELING  OF  THE  BAUSCHINGER  AFFECTED  RANGE 

Note  that  /?is  reported  in  earlier  experimental  work  (ref  7)  relating  to  the  determination  of 
BEF.  The  regime  D-E-F,  interpolated  from  Reference  7,  has  been  fitted  numerically  by  Parker  et 
al.;  see  equation  (2)  in  Reference  8.  In  the  current  notation  this  equation  becomes 

/Gy  =0.064113  +  0.34816  exp(-12.871  Cpt^*) 

+  [1.1619  +  0.5975  exp(-2.175  ept‘*)] 

-  [0.54959  +  0. 74913  exp( -9.637  ept^ *)]  (2) 

Here  all  strains  are  expressed  as  percentage  strain  and  Cpt^*  is  the  maximum  percentage 
plastic  strain  during  prior  tensile  loading.  gJ^^  and  are  total  compressive  stress  and 
compressive  strain,  respectively,  during  cycle  Ul,  i.e.,  is  measured  horizontally  from  point  D 
in  Figure  1  and  Gc^^  is  measured  vertically.  If  G^^  exceeds  the  elastic  value  defined  by  D-E,  the 
elastic  value  is  assumed. 

Equation  (2)  performs  extremely  well  for  current  gun  steels  of  up  to  1100  MPa  nominal  yield 
strength  and  for  the  ranges  of  unloading  plastic  strain,  ep^\  experienced  during  hydraulic 
autofrettage.  However,  the  equation  is  quadratic  in  form  and  is  therefore  unable  to  model  the 
asymptotic  behavior  at  high  compressive  strain  shown  as  a  hatched  line  in  Figure  1.  The 
nonlinear  kinematic  hardening  model  could,  in  principle,  be  used  to  model  over  E-F  during 
the  unloading  cycle  Ul.  This  model  takes  the  form  (ref  1) 

Gc^^  =  atanhix.  epi^^)  +  fC  (3) 

where  a,  Xj  and  Arare  constants  for  the  particular  material. 


4 


Equation  (3)  does  exhibit  an  asymptotic  behavior  at  high  compressive  strain  levels. 
Unfortunately,  since  X'is  analogous  to  ^Gy,  and  P  is  itself  a  function  of  maximum  initial  plastic 
strain,  ept‘*,  this  form  of  the  nonlinear  kinematic  hardening  model  cannot  incorporate 
Bauschinger  effect.  Equation  (4)  is  proposed  for  improved  modeling  of  the  nonlinear  region  E-F. 
It  incorporates  Bauschinger  effect  and  ensures  correct  asymptotic  behavior,  including  strain 
hardening,  at  high  reverse  strains.  Note  that  both  /and  ^are  assumed  to  be  functions  of 
maximum  initial  plastic  strain,  Cpf*-,  the  form  of  each  of  these  functions  is  reported  in  upcoming 
sections. 


=  {[1  +  a- P].  tanhiy.  epP‘)}  +  P+d.  CpP^  (4) 

Since  equation  (4)  must  encompass  the  possibility  that  the  uniaxial  specimen  is  not  deformed 
plastically  in  tension  prior  to  compressive  loading  (i.e.,  Cpt^*  =  0  giving  p  =  1),  it  is  a 
requirement  that 

jicpt'*  =  0)  =  c  (5) 


which  condition  is  enforced  in  upcoming  numerical  fitting  procedures. 

Rees  (ref  9)  argues  that  the  only  meaningful  definition  of  P  is  limit  of  proportionality,  rather 
than  the  usual  choice  of  offsets  such  as  0.2, 0. 1,  or  0.05%  in  the  case  of  Reference  7.  With  the 
benefit  of  hindsight,  it  is  clear  that  a  very  low  value  of  offset  is  desirable  in  determining  p.  An 
offset  value  of  0.01%  was  adopted  in  analyzing  experimental  data  that  then  corresponds 
reasonably  well  with  the  implicit  values  of  P  for  A723  1100  MPa  steel  determined  via  equation 
(1)  from  its  intercept  with  the  equivalent  elastic  behavior. 

One  potentially  attractive  method  of  determining  /would  be  to  impose  continuity  of  slope 
E^'  at  onset  of  Bauschinger  effect  at  CpP^  =  0,  by  differentiating  equation  (4)  with  respect  to 
,  and  thence  calculating  /directly.  This  option  was  tested  but  rejected  because  it  did  not 
provide  the  required  high  quality  of  fit  over  the  crucial  range  following  onset  of  reversed 
yielding. 

EXPERIMENTAL  RESULTS  AND  NUMERICAL  FITTING 

A723  Steel 

Figure  3  shows  a  typical  experimental  stress-strain  plot  from  a  uniaxial  test  on  A723  1130 
MPa  steel  to  1%  total  initial  strain  (ref  3),  during  LI  and  U1  cycles.  Figure  4  shows  similar  data 
for  3%  total  initial  strain.  The  principal  features  described  in  the  introduction  are  evident. 
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Figure  3.  Experimental  uniaxial  stress-strain  data  compared  with  fitting  to 
equation  (2)  and  equations(l)  and  (4);  1%  total  initial  strain,  A723-1 130  MPa  steel. 


Figure  4.  Experimental  uniaxial  stress-strain  data  compared  with  fitting  to 
equation  (2)  and  equations  (1)  and  (4);  3%  total  initial  strain,  A723-1130  MPa  steel. 
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Figures  3  and  4  also  show  the  fit,  based  upon  equation  (2),  to  the  behavior  during  reversed 
yielding.  Clearly  equation  (2)  models  very  well  the  profile  for  the  strain-ranges  encountered 
during  autofrettage  (Table  1),  but  deviates  at  higher  reversed-strains.  Figures  3  and  4  also  show 
the  fit,  based  upon  equations  (1)  and  (4).  y5(0.01%  offset)  was  calculated  from  experimental  data 
and  then  fitted  as  a  function  of  Cpf*.  y  was  then  obtained  as  a  function  of  ept‘*  from  a  "best  fit" 
to  the  full  set  of  individual  experimental  profiles  for  each  of  the  steels  examined,  incorporating 
the  constraint  imposed  by  equation  (5).  Table  2  shows  the  numerical  fit,  including  values  used  in 
equations  (1)  and  (4),  together  with  mean  values  of  material  properties. 

Table  2.  A723  Steel  Numerically  Fit  to  Experimental  Data 
(Note  that  epi^‘*  is  percentage  value.) 


a  =  0,c  =  0,d  =  0.0l3  A723-1130MPa 

a  =  0,c  =  0,d  =  0.035  A723- 1 380  MPa 

eUI/eLi  =  1 . 0.15  Tanh{1.2  Cpf*) 


=  0. 1 684  [7an;2(  1  -  Cpt'  +  0.11  Cpt'  <  =  1 
P=0.\1  e/‘>l 

r=  1 .2858  {ept' A723- 1 1 30  MPa 
y=  1.1244  A723-1380  MPa 


Material  Properties  (Mean  Values) 

Gy  (0. 1  %)  =  1 139  MPa,  Eu  =  209  GPa  A723-1 1 30  MPa 

Gy  (0.1%)  =  1344  MPa,  Eu  =  203  GPa  A723-1380  MPa 


HY180  Steel 


In  the  case  of  HY180  steel,  the  nonlinearity  of  the  initial  plastic  loading  phase  requires  an 
appropriate  numerical  fit,  whereas  the  linear  strain  hardening  is  zero.  Furthermore,  0.05%  offset 
definition  of  initial  yield  is  required  to  obtain  a  satisfactory  fit.  Table  3  shows  the  numerical  fit, 
including  values  used  in  equations  (1)  and  (4),  together  with  mean  values  of  material  properties. 

Table  3.  HY180  Steel  Numerically  Fit  to  Experimental  Data 

(Note  that  ept^*  is  percentage  value.) 


a  =  0.23\,c  =  2,d  =  0 
^ui/^Li  =  1 . 0.17  Tanh(l.2  ept‘*) 


p=\-  0.832  Tanh{0.%5  Cpt^*) 
r=  2  -  1.3  Tanh{0.6  Cpt'*) 


Material  Properties  (Mean  Values) 

Gy  (0.05%)  =1169  MPa,  E^‘  =  201 .9  GPa 
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Figures  5  and  6  show  experimental  results  for  1 .  and  4%  total  strain  compared  with  proposed 
numerical  fit. 


Figure  5.  Experimental  uniaxial  stress-strain  data  compared  with  fitting  to 
equations  (1)  and  (4)  and  to  modulus  ratio;  1%  total  initial  strain,  HY180  steel. 


Figure  6.  Experimental  uniaxial  stress-strain  data  compared  with  fitting  to 
equations  (1)  and  (4)  and  to  modulus  ratio;  4%  total  initial  strain,  HY180  steel. 


PH  13-8  Mo  Steels 


In  the  case  of  PH  13-8  Mo  steels,  there  is  nonlinearity  of  the  initial  plastic  loading  phase  that 
requires  an  appropriate  numerical  fit,  and  the  0.05%  offset  definition  of  initial  yield  is  again 
required  to  obtain  a  satisfactory  fit.  Strain  softening  is  evident.  Table  4  shows  the  numerical  fit, 
including  values  used  in  equations  (1)  and  (4),  together  with  mean  values  of  material  properties. 

Table  4.  PH  13-8  Mo  Steels  Numerically  Fit  to  Experimental  Data 
(Note  that  epi^^*  is  percentage  value.) 


o  =  0.085,  c  =  3.9,  = -0.009 
_  j  _  0.09221  ept'*  +  0.0083  {ept'*f 


/^  =  1  -  0.85*7an;i(4  Cpt'*) 
y=3.9  -3m*Tanh{l>*  Cpt'*) 


Material  Properties  (Mean  Values) 
fTy  (0.05%)  =  1346  MPa,  =  206  GPa  PH  13-8  Mo 
Cy  (0.05%)  =  1302  MPa,  =  195  GPa  PH  13-8  Mo  Super-Tough 


Figures  7  and  8  show  experimental  results  for  2  and  3%  total  initial  tensile  strain  compared 
with  proposed  numerical  fit.  Note  that  two  types  of  PH  13-8  Mo  steel  were  examined.  While  the 
expressions  in  Table  4  adequately  represent  the  normalized  fit  for  both  types,  there  was  a 
significant  difference  in  the  mean  values  of  0.05%  offset  yield  strength.  This  effect  will  be  of 
significance  in  future  work  involving  calculation  of  autofrettage  residual  stresses  for  the  two 
types  of  PH  13-8  Mo  steel. 
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Figure  7.  Experimental  uniaxial  stress-strain  data  compared  with  fitting  to 
equations  (1)  and  (4)  and  to  modulus  ratio;  2%  total  initial  strain,  PH  13-8  Mo  steel. 


Strain  (%) 


Figure  8.  Experimental  uniaxial  stress-strain  data  compared  with  fitting  to  equations  (1) 
and  (4)  and  to  modulus  ratio;  3%  total  initial  strain,  PH  13-8  Mo  super-tough  steel. 
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CONCLUSIONS  AND  FUTURE  WORK 


Gun  steels  have  been  subjected  to  a  series  of  uniaxial  tests.  These  tests  are  intended  to 
provide  equivalent  stress  data  for  input  to  subsequent  numerical  analyses.  The  steels  are  A723- 
1130  MPa  (nominal)  yield  strength  and  possible  alternatives,  namely  A723-1380  MPa  (nominal) 
yield  strength,  HY80-1180  MPa  (nominal)  yield  strength,  PH  13-8  Mo  and  PH  13-8  Mo  super¬ 
tough. 

For  each  steel  the  experimental  data  during  initial  and  reversed  loading  can  be  fitted  with 
considerable  accuracy  using  a  new  variant  of  the  nonlinear  kinematic  hardening  model  that 
includes  Bauschinger  effect  as  a  function  of  prior  tensile  plastic  strain.  This  new  variant  provides 
an  excellent  fit  within  the  strain  ranges  experienced  in  real  gun  tubes  during  autofrettage.  It  also 
provides  the  desired  asymptotic  behavior  at  much  higher  reverse-strain  levels  and  therefore 
represents  an  improvement  over  previous  models. 

Future  work  will  encompass: 

•  Calculation  of  residual  stress  profiles  and  associated  fatigue  lifetimes  in  tubes  made  from 
the  candidate  steels. 

•  Prediction  of  onset  of  reyielding  during  proof  testing. 

•  Detailed  prediction  of  nonlinear  strain  behavior  during  typical  proof  testing  and  firing. 
This  may,  in  turn,  lead  to  a  reassessment  of  the  safe  maximum  pressure  calculation  for 
autofrettaged  tubes. 
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NOMENCLATURE 


a,  c,  d 
^  UJ 


epi 


Ll 


^pl 


LJ* 


^pl 

eLI 

j?Vl 


UI 


P 

r 


Cpi 


Ll 


.UI 


(^pl 


UI 


(Ty 


Coefficients  in  equation  (1)  relating  to  initial  plastic  deformation 
Total  compressive  strain  during  reversed  loading 
Plastic  strain  during  initial  plastic  deformation 
Maximum  plastic  strain  during  initial  plastic  deformation 
Nonlinear  strain  during  reversed  loading 
Elastic  modulus  during  initial  tensile  loading 
Elastic  modulus  during  initial  load  reversal 
Traditional  nonlinear  kinematic  hardening  coefficients 
Bauschinger  effect  factor,  function  of 

I  I  ^ 

Proposed  nonlinear  kinematic  hardening,  function  of  Cpi 
Stress  during  initial  plastic  deformation 
Total  compressive  stress  during  reversed  loading 
Nonlinear  stress  during  reversed  loading 
Material  yield  stress  obtained  from  given  offset 
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